Evolution of Magnetic Fields in Freely Decaying Magnetohydrodynamic Turbulence 
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We study the evolution of magnetic fields in freely decaying magnetohydrodynamic turbulence. By 
quasi-linearizing the Navier-Stokes equation, we solve analytically the induction equation in quasi- 
normal approximation. We find that, if the magnetic field is not helical, the magnetic energy and 
correlation length evolve in time respectively as Eb oc t^^^^+^'/t^+p) ^ ^2/{3+p) ^ .^^^here p is the 

index of initial power-law spectrum. In the helical case, the magnetic helicity is an almost conserved 
quantity and forces the magnetic energy and correlation length to scale as Eb oc {logt)^^^t~'^^^ and 

^B oc {\0gt)-'/h'^^ 
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The origin of presently-observed large scale magnetic 
fields throughout the universe is still unclear Essen- 
tially, there are two possible classes of mechanisms to 
produce cosmic fields depending on when they are gen- 
erated: Astrophysical mechanisms acting during or after 
large-scale structure formation, and mechanisms acting 
in the primordial universe. Magnetic fields created in the 
early universe (except those generated during infiation), 
usually suffer from a "small-scale problem" , that is their 
comoving correlation length is much smaller then the 
characteristic scale of the observed cosmic fields. How- 
ever, if magnetohydrodynamic (MHD) turbulence oper- 
ates during their evolution, an enhancement of correla- 
tion length can occur, especially if the magnetic field is 
helical. As pointed out by Banerjee and Jedamzik 
the evolution of a magnetic field in the early universe 
goes through different phases depending on the particu- 
lar conditions of the primordial plasma. In this paper, 
we are interested in the case of magnetic fields evolving 
in the turbulent primordial universe well before recombi- 
nation epoch and when kinematic dissipative effects are 
due to diffusing particles. Therefore, we are concerned 
with the so-called phase of "turbulent MHD". In other 
phases, such as "viscous MHD" and "MHD with ambipo- 
lar diffusion" described in Ref. [1], the dynamics of the 
magnetic field is very different from that studied here. 
The problem of determining the evolution properties of 
magnetic fields in MHD turbulence has been deeply and 
widely discussed in the literature using different methods 
and approximations. A direct integration of the full set 
of MHD equations would allow us to deeply understand 
the dynamics of freely decaying MHD turbulence. How- 
ever, MHD equations are quite difficult to handle due to 
their high non-linearity and it has not been yet brought 
in a definitive verdict for the evolution laws of magnetic 
energy and correlation length (for recent numerical stud- 
ies of freely decaying magnetohydrodynamic turbulence 
see, e.g., Ref. 0, H, 

The turbulent MHD equations for incompressible flu- 
ids, in the case of non-expanding universe, are 



i9tv -I- V Vv + Vp - i/V^v = J X B, (1) 
= V X (v X B) + ryV^B, (2) 

and V • V = V • B = 0. Here, v is the velocity of bulk fluid 
motion, B the magnetic fleld, J = V x B the magnetic 
current, v the kinematic viscosity, r/ the resistivity. The 
thermal pressure of the fluid, p, is not an independent 
variable since, taking the divergence of Eq. (jlj, it can be 
expressed as a function of B and v. 
In the case of expanding universe in the radiation era, it 
has been shown that the MHD equations are the same 
as Eqs. (II])-(I2]) provided that time, coordinates, and dy- 
namical variables are replaced by the following quanti- 
ties (see, e.g., Ref. [2]): t ^ i ~ Ja^^dt, x — > x = ax, 
B — !■ B = a^B, v ^ 9 = a~^v, r] fj = a~^r], where 
a{t) is the expansion parameter. Because of the formal 
coincidence of the MHD equations in the expanding and 
non-expanding universe, we can study the evolution of 
magnetic fields in MHD turbulence in both cases in a 
similar way. For definiteness, in this paper we shall con- 
sider only the case of non-expanding universe. 
It is useful to define the kinetic and magnetic Reynolds 
numbers. Re = vl/v and Re^ = vI/Vj where v and 
I are the typical velocity and length scale of the fluid 
motion. Magnetohydrodynamic turbulence occurs when 
Re 3> 1 and Re^ 3> 1. We are interested in the evolu- 
tion of statistically homogeneous and isotropic magnetic 
fields. This means that the two-point correlation tensor 
Cij(x,y) = {Bi{x)Bj{y)), where (...) denotes ensemble 
average, is a function of |x — y| only and transforms as 
an SO{3) tensor. In terms of the Fourier amplitudes of 



the magnetic field, i?i(k, t 
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{B,{k,t)B,{p,t)) = [(27r)72]5(k + p) 

X (Sij - kikj)S{k,t) + i£ijkkkA{k,t) , (3) 

where ki — ki/k, k — |k|, and Sijk is the totally antisym- 
metric tensor. The functions S and A denote the sym- 
metric and antisymmetric parts of the correlator. They 
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are related to the magnetic energy and helicity densities 
in the volume V through EB{t) = (1/2F) Jy d^x (B^) = 
J^dk£Bik,t) and H Bit) = il/V)Jyd^x{A-B) = 
J^dknBik,t), where £b = k'^S/{2TT)^ and Hb = 
kA/{27T'^) are the magnetic energy and magnetic helic- 
ity density spectra, and A is the vector potential. The 
kinetic energy, Ey (t) , is defined as the magnetic one with 
B replaced by v. We remember that for all magnetic 
field configurations, the magnetic helicity spectrum must 
satisfy the "realizability condition" Q: |H_b| < 2k^^£B- 
The magnetic field is said to be "maximally helical" if, 
for all fc, TLb is of the same sign and saturates the above 
inequality. Moreover, the magnetic helicity is conserved 
when 77 = since Q 9*^3 = -{2ri/V)jyd^x {3 ■ B) = 
—2riJ^dkk^TiB- The relevant length scale in MHD the- 
ory, the so-called magnetic correlation length, is the char- 
acteristic length associated with the large magnetic en- 
ergy eddies of turbulence and is defined by: £,B{t) ~ 
Eg^ dkk~^ £b- The integral form of the realizabil- 
ity condition takes the form: \Hb\ < 2£^bEb- 
Since we are interested in the case of large Reynolds num- 
bers, we neglect the dissipation term in Eq. More- 
over, as in Ref. f^, we quasi-linearize the Navier-Stokes 
equation ([1]) neglecting the quadratic term (v • V)v. ^ 
This corresponds to neglect small scale components of 
velocity field and to assume that the Lorentz force, 
Fl = J X B, acting on the charged particles of the fluid 
"drives" the development of turbulence on larger scales: 
9tV ~ Fi. Although the validity of this approximation 
can be verified only by a numerical analysis, its use is 
justified a posteriori since our results, as we will find, are 
in agreement with a numerical simulations of full MHD 
equations performed in Ref. f3|. Finally, we make the 
common approximation, 

v~TdFi, (4) 

where the "drag time" is the fluid-response time to the 
Lorentz force introduced by Sigl in Ref. 8] . We note that 
Td and the characteristic time associated with kinetic tur- 
bulence, the so-called eddy turnover time Teddy — l/v, are 
related by ~ Freddy, since 1 ~ |v|/|TdFL| ~ TTeddy/rd- 
Here, F is the ratio of the kinetic and magnetic energy. 
What is observed in numerical simulation of MHD equa- 
tions is that, in the non-helical case, turbulence pro- 
ceed toward a state of equipartition between magnetic 
and kinetic energies (F ~ 1) while, in the helical case, 
though there is no evidence of equipartition, the ratio F 



^ If we decompose, in the spirit of mean- field-theory the ve- 
locity field into an (almost uniform) average part and a weak, 
small-scale fluctuating part, v = v-|-<5v with |<5v| <g |v|, we have 
V • Vv ~ V • V5v. Comparing this term with the Lorentz force, 
we get |v • Vv|/|Fl| ~ r|(5v|/|v|, where T = E^/Eb- Hence, 
the quasi-linear approximation is valid as long as the condition 
r 3> |(5v|/|v| is satisfied. 



approaches asymptotically to a constant value. 
Inserting the above expression for v into the induction 
equation we get in Fourier space: 

/d^ V f d^ Q 
l2nf J (2^ eyfefcj9rS,(q) x 

[skrsBn (p- q)B„ (k- p) -e^sm-Bfe {k-p)Bm (p- q)] , (5) 

where summation over repeated indexes is understood. 
We will work in "quasi-normal approximation" and sup- 
pose that the four-point correlator can be decomposed, 
in terms of two-point correlator, as Q: 

{B,{k)Bj{p)Bk{q)Biir)) 

= (i?.(k)i?,(p))(Bfe(q)i?,(r)) 

+ {B,{k)Bk{q)){B,{p)Bi{r)) 

+ {Bi{k)Bi{r)){B^{p)Bk{ci)). (6) 

Multiplying Eq. ([5]) respectively by i?i(k) and Ai(k), and 
then averaging out we arrive at the following equations 
for the magnetic energy and helicity spectra: 

dt£B = -2r/efffc2£s + aBk'^HB, (7) 
dtHB = -2?7efffc^7is + iaBSs, (8) 

where we have introduced ?7off(^) = ^7 + '^^EBTd/^i and 
ctsit) = —HBTd/i'ir]). For simplicity, we will restrict our 
analysis to magnetic fields with initial "fractional helic- 
ity": HBik,0) = /iB7i:g^^(fc,0), where < /is < 1 is 
the fraction of the initial maximal helicity Ti-B^^ikji) — 
2k~^£B{k,t). In this case, the solution of Eqs. (I?!)-® is: 

£B{k,t) ^ £B{k,0)exp{-2k^eU 

X [cosh(2Ma) -I- hB sinh(2Ha)], (9) 
HBik, t) = HT^'ik, 0) exp{-2kHlJ 

X [sinh(2Ha) + Hb cosh(2M„)], (10) 

where we have defined the "dissipation" and "alpha" 
lengths, -^disslO = Jodtrjcs and ia(t) = J^dtaB- From 
Eqs. ([Sjl- fTU]) we immediately get that magnetic fields 
with maximal initial helicity, /i^ = 1, remain maxi- 
mally helical for all times: Hb — 2k^^£B- To pro- 
ceed further, we assume that the initial magnetic en- 
ergy spectrum can be represented by the following sim- 
ple function: £B{k,0) = Asfc^ exp(— 2fc^^^), where 
and Ib are constants. For k <C the magnetic en- 
ergy spectrum possesses a power law behavior, while 
for large k it is suppressed exponentially in order to 
have finite energy. The exponential cut-off, £b, is re- 
lated to the initial correlation length by £b — £,b{0)/Cb, 
where (b = V2T{p/2)/T[{l+p)/2] and T{x) is the Euler 
gamma function. In Ref. it was shown that ana- 
lyticity of the correlator (x,y) defined on a compact 
support forces the spectral index p to be even and equal 
or larger than 4. Now, inserting Eqs. (f9|l- (fT0|) into the 
expressions for the magnetic energy and helicity we find: 
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where iFi{a,b; z) is the Kummer confluent hypergeo- 
metric function, and we have defined Cdiss = ^diss/^B, 
Ca = C/^B, and X = Ca/(1 + Cdiss)'^'- Equations 
(fTT|) - (|12p are integral equations for the magnetic en- 
ergy and helicity. They can be solved once the ex- 
plicit expression for the drag time is given. This can 
be done if we consider the scahng properties of the in- 
duction equation. It is well-known that the full MHD 
equations (neglecting dissipative terms) are invariant un- 
der the scaling transformations x ^x, t i^~'^t, 
V ^ ^'"v, B ^ ^''B, where ^ > is the "scaling fac- 
tor" and r is an arbitrary real parameter [lo| . Now, 
imposing that also the "reduced" MHD equations ([2]) 
and (III) are invariant under these scaling transforma- 
tions, we get that the drag time is linear in time. Tak- 
ing into account the relation between Td and Teddy pre- 
viously discussed, we also have that the eddy turnover 
time is asymptotically linear in time. This allow us to 
write the drag time as Td{t) ~ r(0)[Tcddy (0) + 7^], where 
7 = [r(oo)/r(0)] limt_>oo Teddy (i)/i IS a constant, whose 
explicit value is inessential for the following discussion. 
It is useful to define accurately the magnetic Reynolds 
number and the eddy turnover time: Re^ = Wims^s/jy, 
and Teddy = Cs/wrms, whcrc, as typical length scale 
and velocity, we used the magnetic correlation length 
and the root-mean-square value of the velocity field, 
= {l/V)Jyd^x{-v'^) = 2Ey. With the aid of the 
above definitions and introducing the normalized time 
T = t/Teddy(0), thc integral equations (fTTj) and (fT2)) can 
be transformed into the differential equations 



dr 



ReB{0) 3 



+ oC|(1 + 7t) 



Eb{t) 
Eb{0) 



(13) 



dT=4^-'^-^^-(°)(^+^^)^i$y'(^^) 

where Eb and Hb, as a function of Cdiss and Cq., are 
given by Eqs. (fTTj) and (fT2|) . For large magnetic Reynolds 
numbers, the first term in the left-hand-side of Eq. (|13p 
can be neglected with respect to the second one. 

In the non- helical case, Hb = 0, the solution of 
Eqs. ((I3|)-(lll|) is Ca = 0, that is HB{t) = for aU 
times, and CLs = [1 + Kdiss(2T + r2)]2/(3+P) _ ^here 
Kdiss = 7(3 + p)Cb/6- This, in turn, gives for t » 1: 

Eb{t) ~ keEb{0) T-2(i+P)/(3+f), (15) 
^B{r) ^ n^UOy^^'^'\ (16) 
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FIG. 1: Magnetic energy spectrum in the non-helical case 
for p = 4, with 7=1. The dotted line corresponds to the 
initial spectrum, while continuous lines correspond, from left 



to right, to t/Teddy(0) = 1, 10, 10^ 



,10' 



where .b = ^dil:^^^^^''"^^ -nd = ^^^^^ . 
It is interesting to observe that, starting from self- 
similarity of MHD equations, Olesen obtained the fol- 
lowing expression for the magnetic energy spectrum [id ] : 

2 

EB{k,t) — XBkP ipBikt^+p), where Xb is a constant, ■0b 
is an unknown scaling-invariant function, and p is the 
power-law exponent of the initial magnetic energy spec- 
trum. Our approach to MHD equations fixes the ex- 
pression of the scaling-invariant function to ^b{x) = 
exp[-2(a;/x,)2], with = K^-^£BMTeddy(0)]2/(3+P). 
In Fig. 1, we plot the spectrum of the magnetic energy for 
the case p = 4 at different times. It is clear that, in the 
non-helical case, the decay of the magnetic field and the 
growth of the correlation length proceed through the so- 
called selective decay discussed by Son in Rcf. lli] : there 
is no direct transfer of magnetic energy from small scales 
(large wavenumbers) to large scales (small wavenumbers) 
but, simply, modes with larger wavenumbers decay faster 
than those whose wavenumbers are small. Consequently, 
as the turbulence operates, the magnetic field survives 
only on larger and larger scales. 

In the helical case, the evolution of the system goes 
through two different regimes depending on the value 
of x{t) which is an increasing function of time. Un- 
til when X ^ 1 the system behaves as if the magnetic 
helicity were zero: the system evolves by selective de- 
cay and, consequently, the asymptotic solutions are the 
same obtained previously. Afterwards, when x S> 1, 
the system enters and persists in a phase characterized 
by a transfer of magnetic energy from small to large 
scales, a mechanism known as inverse cascade [5]. The 
asymptotic (t 00) solutions in this latter phase are: 
Cdiss(T) ~ Cdiss(lnT)i/6T2/3 and Ur) ~ c^ln T)2/3r2/3, 

where Cdis. = (7 V12)'/'Cb/^s^' and c„ = {Ap/Sy/^diss- 
Consequently, we have: 

EBiT) ~ ceEb{0) (lnT)i/3 T-2/3, (17) 

CB(T)c.c«eB(0)(lnT)-i/V2/3, (18) 

with ce = (2^/37)^/^/1^'''^, and = hsc]^^. From 
the above equations, we directly obtain the relation 
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FIG. 2: Result of a numerical integration of Eqs. (|13p - p4p for 
Res(O) = 10^^, p = 4, ft = 10"^, with 7 = 1. Upper panel: 
magnetic energy; middle panel: correlation length; Dotted 
lines correspond to analytical expansions. Lower panel: mag- 
netic energy spectrum; the dotted line corresponds to the ini- 
tial spectrum, while continuous lines correspond, from left to 
right, to t/roddy(0) = 1, 10, 10^ 10^ 



Eb^b — Hb/2. This means that a magnetic field with 
initial fractional helicity becomes maximally helical ap- 
proximatively after the system enters into the inverse- 
cascade regime. More accurately, we can find the time 
when this happens, t^, matching the product of asymp- 
totic solutions (Unil-ITCl) and ((T71)-(ITH]). It resuhs: Th ~ 
'^diss^ hj^'^^^^^'^^ ■ In Fig. 2 we present the result of a nu- 
merical integration of Eqs. (Unil-lllll) for ReB(O) = 10^^, 
p = 4, and h = 10"'^. It is evident from the figure that 
the analytical expansions [non-helical solution for t 
and Eqs. P?)l -(fT5 |l for t > Th] fit very well the numeri- 
cal solution. Because of quasi-conservation of magnetic 
helicity, small-scale modes are not dissipated during the 
decay but their energy is transferred to larger scales: this 
process of inverse cascade is manifest in the magnetic en- 
ergy spectrum shown in Fig. 2. 

It is worth noting that our final results, Eqs. P^ - P^ . 
apart logarithmic factors, agree very well both with scal- 
ing arguments and results of a numerical integration of 
full MHD equations presented in Ref. 0. 

In conclusion, we have studied the evolution of statis- 
tically homogeneous and isotropic magnetic fields in the 



context of freely decaying magnetohydrodynamic turbu- 
lence. By quasi-linearizing the Navier-Stokes equation, 
we have solved analytically the induction equation in 
quasi-normal approximation. We have found that, if the 
initial magnetic field is not helical, the evolution of the 
magnetic field proceeds through selective decay of mag- 
netic modes: magnetic power on small scales is washed 
out by turbulence effects more effectively than on large 
scales. During this process, the correlation length grows 
as oc i^/(3+p)^ while the magnetic energy decays in 
time as Eb oc t~^(^+P)/(^+P\ where p is the index of the 
initial power-law spectrum. In the helical case, the evo- 
lution of the system goes through two different phases: 
selective-decay phase in which the system evolves as if the 
magnetic helicity were zero and inverse-cascade phase. 
The first phase ends when quasi-conservation of mag- 
netic helicity starts to trigger an inverse cascade of the 
magnetic field: small-scale modes are no more completely 
dissipated during turbulence but their energy is partially 
transferred to larger scales. This causes a faster growth 
of the correlation length and a slower dissipation of the 
magnetic energy with respect to the non-helical case. 
The time when the system enters into the inverse-cascade 
regime is proportional to hg^^~^^^^^^ times the initial eddy 
turnover time, where Hb is fraction of the maximal initial 
magnetic helicity. Moreover, the process of inverse cas- 
cade erases any information about the initial structure 
of the magnetic field, so that the evolution laws of en- 
ergy and correlation length are Eb oc (logi)^/'^^"^/'^ and 
^b oc (logi)~^/'^t^/'', whatever is the value of p. 
In a cosmological context, these results are of interest 
when studying the evolution of primordial magnetic fields 
before neutrino decoupling. Indeed, during the period 
of neutrino (or photon) free-streaming, as well as after 
recombination, the equations governing the evolution of 
magnetic fields differ from those studied here and then 
our results do not apply. Nevertheless, our approach to 
MHD equations can be suitably extended to these last 
cases and an appropriate analysis is in progress. 
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